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We show that in the metalhc phase of a two dimensional electron gas the spin-orbit coupling due 
to structure inversion asymmetry leads to a characteristic anisotropy in the magnetoconductance. 
Within the assumption that the metallic phase can be described by a Fermi liquid, we compute 
the conductivity in the presence of an in-plane magnetic field. Both the spin-orbit coupling and 
the Zeeman coupling with the magnetic field give rise to two spin subbands, in terms of which 
most of the transport properties can be discussed. The strongest conductivity anisotropy occurs for 
Zeeman energies of the order of the Fermi energy corresponding to the depopulation of the upper 
spin subband. The energy scale associated with the spin-orbit coupling controls the strength of the 
effect. More in particular, we find that the detailed behavior and the sign of the anisotropy depends 
on the underlying scattering mechanism. Assuming small angle scattering to be the dominant 
scattering mechanism our results agree with recent measurement on Si-MOSFET's in the vicinity 
of the metal-insulator transition. 



I. INTRODUCTION 

In recent years the transport properties of two- 
dimensional (2D) electron systems in Si-MOSFETs and 
semiconductor heterostructures have been the subject of 
a great deal of theoretical and experimental activity. One 
reason of this is the possibility of a metal-insulator tran- 
sition (MIT) as the density of the 2D system is varied 
101 . A "critical" density identified by a weakly temper- 
ature dependent resistivity separates an "insulating" re- 
gion from a "metallic" one. In the latter, upon lowering 
the temperature roughly below IK, the resistivity drops 
by almost one order of magnitude, whereas in the former 
the resistivity increases quickly. The origin of the MIT is 
still unclear and there is a wide debate in the literature 
concerning the relevance of the possible mechanisms. For 
recent reviews on the subject one may see Rcf. and 

i- 

On the insulating side the resistance resembles Efros- 
Shklovskii hopping which points towards the rele- 
vance of localized states and a soft gap in the density 
of states; also the thermopower has been interpreted in 
terms of an Anderson insulator From this point of 
view Anderson localization seems to be relevant for the 
transition. Whereas the standard single parameter scal- 
ing theory predicts no Anderson transition in two dimen- 
sions, on the basis of the existing scaling theory of inter- 
acting electrons a transition cannot be excluded 
Besides, several properties of the metallic side, including 
the low density, i.e., the large value of the electron gas 
parameter Ts , the strong increase of the electron factor 
1^,^, the anomalous magnetoresistance and a 

vanishing compressibility at the critical density [ p^ can 
indeed be seen as hints for the relevance of the electron 
correlations; quantum corrections to the conductivity on 
the other hand seem to be only small ]T6|-p^, at least 
in the range of temperatures where the main drop of the 



resistivity occurs ||T^ . 

The main and still open questions are: (1) Is there a 
real quantum phase transition, or is just a temperature 
dependent crossover phenomenon observed? (2) What is 
the nature of the insulating state? (3) Can the metallic 
phase be described in terms of the standard Fermi-liquid 
theory? 

In this paper we will not touch point one and two. 
Concerning point three, there is accumulating evidence 
that the answer is "yes" : at the lowest temperatures and 
weak magnetic fields standard weak localization physics 
has been observed in a number of different systems. 

In order to obtain a better understanding of the mate- 
rials, different aspects have to be investigated. We study 
the influence of spin-orbit coupling on the magnetocon- 
ductance. 

We assume that the metallic phase can be described 
in terms of the standard Fermi-liquid theory, and there- 
fore we calculate transport properties in the framework of 
the Drude-Boltzmann theory. In Si-MOSFET's with low 
electron density the major spin-orbit term is believed to 
be due to the lack of inversion symmetry of the confining 
potential. Originally, there has been the suggestion that 
the spin-orbit coupling may be responsible for the ob- 
served metallic behavior of the resistivity . Although 
at present such a view is not confirmed by the experi- 
ments, it is certainly useful to assess the relevance of the 
spin-orbit coupling in these systems as well as to under- 
stand whether it is related to the puzzling magnetoresis- 
tance in a parallel field. In this respect, Chen et al ]2l| ] 
pointed out that such a spin-orbit coupling would induce 
an anisotropy in the magnetoresistance. The resistivi- 
ties measured parallel and perpendicular to the in-plane 
magnetic field are different. Such an anisotropy has been 
observed in Si-MOSFETs |2|. 

However, while the original theoretical proposal in Ref. 
pl[ has been formulated in the localized limit, the exper- 
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imental results have been obtained near and across the 
MIT. For this reason a direct comparison between the- 
ory and experiment is difficult. In a recent paper p^ ] 
we have provided a theory valid in the metallic regime. 
We were able to reproduce correctly a number of the 
experimental findings. The only drawback of our the- 
ory was the wrong sign of the effect, and we speculated 
that this could be due to the oversimplified assumption 
of a pure s-wave impurity scattering. It is also worth 
mentioning that anisotropic magnetoresistance has been 
observed in GaAs electron systems [Q, and GaAs hole 
systems [ p5p^ . In these systems, the interpretation is 
complicated by the intrinsic anisotropy associated with 
the crystal structure. 

In this paper, we study how the form of the microscopic 
scattering potential affects the anisotropy of the magne- 
toconductance in the presence of the spin-orbit coupling. 
In particular, we extend our previous analysis to include 
the case of small angle scattering. The layout of our pa- 
per is the following. In section II we introduce the model 
and make a first analysis by means of the Boltzmann 
equation. We adopt the relaxation-time approximation 
assuming a relaxation rate independent from the posi- 
tion on the Fermi surface. In section III, starting from 
the Kubo formula we calculate the conductivity within 
a Green function approach. As a first step, we avoid 
to specify a microscopic scattering mechanism and we 
simply introduce, by hand, a lifetime in the one-particle 
Green function. In so doing we neglect vertex corrections 
but take into account level broadening and inter-band 
transitions. As it will be clear in the following, while 
vertex corrections may give rise to significant changes in 
the conductivity, the inclusion of level broadening affects 
only weakly the results. A more sophisticated analysis is 
carried out in sections IV and V, where we will specify 
a microscopic mechanism responsible for the finite con- 
ductivity, namely elastic scattering. In section IV we will 
assume s-wave scattering. The Green functions and ver- 
tex corrections will be calculated in the self-consistent 
Born approximation. In section V we will explore the 
consequences of a strongly angle dependent scattering po- 
tential. In the limit of strong forward scattering we will 
find that the sign of the anisotropy is consistent with the 
experiments. Finally, section VI will contain our conclu- 
sions. 

II. THE MODEL AND THE RELAXATION TIME 
APPROXIMATION 

We start from the model Hamiltonian 

1 

H^- h acr • p A - -.g^Bcr ■ B (1) 

2m 2 

where a is a parameter describing the spin-orbit coupling 
due to the confinement field a is the Pauli matrices 



vector, and is a unit vector perpendicular to the 2D 
system. As pointed out recently by Winkler [|8|, this 
so-called Rashba model applies for 2D electron systems 
(like Si-MOS in Ref. [|2|), but not for heavy-hole states 
which are relevant e.g. in Si/SiGe quantum wells. To 
simplify the notation we introduce the Zeeman energy 
UJs — \gi^BB (note that this definition differs by a fac- 
tor of 2 from the standard one). The ^-factor in nSi- 
MOS structures is near two but is growing with decreas- 
ing electron density [ pl| . In the presence of spin-orbit 
coupling and magnetic field, the original spin degenerate 
band splits into two bands with energy dispersion 

i^±(p)=pV2m±r!(p) (2) 

f2(p) = ^J{apy - UsY + {apxY- (3) 

We have chosen the direction of the in-plane magnetic 
field as the x-axis. According to the Boltzmann trans- 
port theory within the relaxation time approximation, 
the conductivity tensor at zero temperature is 

where s = ± is the band index and the velocities are 
v± = Vpi?±(p). By further assuming a constant relax- 
ation rate Tp ^ = rtr, the conductivity is proportional to 
the product of the density of states, Ns{fi), times a Fermi 
surface average of the Fermi velocities, 

a,, = eVtr ^ Ns{n){vlvi)FS^. (5) 

s 

We note also that, after integrating by parts Eq.(Q), the 
conductivity may be expressed in terms of the effective 
mass tensor, 

In the absence of spin-orbit scattering this relation be- 
comes simply a = 2e'^T^^n/m, with n the electron den- 
sity. We assume that the density does not change with 
the magnetic field and therefore we have to adjust the 
chemical potential. As long as both bands are occupied, 
the chemical potential is nearly constant, fiiujs) ~ £f- 
For high fields when the upper band is depopulated it 
decreases as /i(ws) ~ 2eF — oJs- The chemical potential 
in the absence of magnetic field is here denoted by Ci;^. 

From Eqs.(|^) and (|^) it is apparent that the trans- 
port time drops when we consider the conductivity rel- 
ative to its value in the absence of magnetic field and 
spin orbit scattering, ctq = c^Nqv^pT^^^. For the same rea- 
son, the relative anisotropy in the conductivity. Act = 
{(^xx — ^yy)l'^Oi does not depend on the scattering time 
Ttr and may be determined from the simple knowledge 
of the band structure. This observation makes clear the 
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main goal of our analysis, i.e. to what extent we can un- 
derstand the small anisotropy in the conductance with- 
out making any assumption about the origin of the strong 
magnetic field dependence of the conductivity itself. 
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FIG. 1. Aa versus magnetic field {ujs — q^bB /2) at fixed 
spin-orbit coupling app = QAep as obtained within the re- 
laxation time approximation with a single scattering rate 

Ttr- 

The evaluation of the conductivity via Eq. (||) or (|^) is 
analytically nasty due to the presence of the square root 
in the energy band dispersion (cf. Eq.(||)). We begin 
then by presenting the result of the numerical integra- 
tion of Eq. (^) in Figure |^ which shows Act for a special 
value of the spin-orbit term as a function of the magnetic 
field. We have also plotted the contributions for the two 
bands separately, and it is seen that at weak magnetic 
field there is a strong cancellation from the contributions 
of both bands. At large field the upper band depopulates 
and therefore does not contribute to the conductivity or 
its anisotropy. 

Some more insight is obtained by considering Eq. (^) in 
the weak and strong field limits where analytical progress 
is possible. First we note that Act must be zero in the 
limits of very weak and very strong magnetic field, when 
the Fermi surface becomes rotational symmetric. When 
ijJs > apF the dispersion of the two bands is approxi- 
mately given by 



E±ip) = _p^/2to ± [ujs - apy + {ap^)'^/2ujs 



(7) 



so that the effective mass in the x-direction, i.e. parallel 
to the magnetic field, is modified according to 



(8) 



For tUs > ef, when only one band is occupied, the 
anisotropy in the conductivity is thus proportional to 
l/nix ~ l/m and explicitly given by 



1 {aPF? 

2 LOsep 



(9) 



For LOg < ep, when both bands contribute, there is a 
partial cancellation: 



^xx ^yy 



1 {appY ^- ~ 

2 ujgep n 



(10) 



1 {c^PP? 
'2 4 ' 



(11) 



since the electron densities in the lower and upper band 
are n± « n(l =p L0alep)/2. 

In order to give the results for the low field limit, we 
first write (0) explicitly for the given band dispersion 



^xx ^yy — C Tt 



{apxY - {^s - aPyY 

f]3 



x[e(M-i?-(p))-e(/.-i?+(p))]. (12) 

When ^ <^ ep one may then approximate the difference 
of the step functions by a delta function which leads to 



(^XX ~ (Tyy _l {appf 



2 6| 

2lT 



(13) 



A9 [app cos{9)]^ ~ [uog — app sin(0)]^ 



27r 02 
Evaluating the integral one arrives at 
CTxx - a-yy 1 (app)'^ 1 - (ujs/app)'^ 



CTO 



o 2 ^ / ^2 ^i^Jo^PF - 1) 
2 ej, [uJslappY 



(14) 



It is interesting to note that ( [T^ ) is identical to the Fermi 
surface average of the anomalous part of the velocity 
v± = p/m ± aep, 

CTxx - CTyy = 2e^ iVoTtr ( (Gp ■ Bx)'^ - (Cp • eyf)ps (15) 

with Bp = (ap — uJsey)/Q.. At low magnetic field the 
vector Bp makes a full rotation, when averaging over the 
Fermi surface so that the two dot products averaged over 
the Fermi surface, in Eq.(^5|), cancel each other. The 
anisotropy of the conductivity vanishes. At large field Bp 
becomes locked in y-direction. The conductivity becomes 
then anisotropic at a magnetic field value uJs = app. 
Fig.^ shows the numerically determined anisotropy for 
different strengths of the spin-orbit coupling. The edge 
behavior SlI app — lOs predicted by Eq.(^4|) is clearly 
visible. 
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FIG. 2. Act versus magnetic field as obtained within the 
relaxation time approximation for various strengths of the 
spin-orbit coupling. 
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When tp > ^ apF the anisotropy is nearly con- 
stant, which agrees both with our weak and strong field 
expansion. 



III. THE CONDUCTIVITY WITHIN THE KUBO 
FORMULA APPROACH 

In this section we go beyond the Boltzmann equation 
approach in the relaxation-time approximation and cal- 
culate the conductivity via the Kubo formula. The start- 
ing point is the response kernel Qij which at zero tem- 
perature is given by 

/°° df 
—Tr[f{p,q)G{p+)J\p,q)G{p^)] 



(16) 



In the above we are using four vectors p = (e, p) and 
q = (w, q) with p± = p± q/2; the trace is over the spin. 
In the presence of the spin-orbit coupling the bare cur- 
rent vertex differs from the standard expression and is 
given by 



J = aa Ae^ 

ni 



(17) 



The evaluation of the dressed vertex will be considered in 
section IV. The conductivity tensor is obtained by con- 
sidering the limit 



~ lim 



(18) 



By splitting the energy integration into regions where the 
Green functions have defined analytical properties, one 
may separate the kernel into two parts. In the first both 
the Green functions are retarded or advanced. This part 
cancels exactly the diamagnetic term, i.e., 

/°° df N 
-Tr [j\p)G{p).P[p)G{p)] - -5,, = 0. (19) 

The cancellation is made explicit in appendix ^ by ex- 
ploiting the Ward identity due to the gauge invariance. 
The second term is given by 

= [23\p)G^ip)JU^)G^ip) 



f{p)G^{p)Ji,^{p)G''{p) 
f{p)G^ip)JUp)G^ip) 



(20) 



The dressed vertex depends whether it is connected to a 
pair of retarded and advanced Green functions or a pair 
of Green functions with equal analytic properties. 



In order to make contact with the analysis of the pre- 
vious section, we temporarily make no assumption re- 
garding a specific scattering mechanism and simply take 
a (retarded) Green function of the form 



(Gf)-i(p,e)=e-£;±(p)+M + i/2T, 



(21) 



where the scattering time t is of some unknown origin. 
Let us also neglect the vertex corrections, i.e. J j- We 
are aware that neglecting vertex corrections often leads 
to serious errors in the calculation of the conductivity, 
but their calculation requires a specific scattering mech- 
anism. The conductivity is then given by 



^■'^■ = -E E («lf(p)l&)Ini(Gf(p))(6|f (p)|a)Im(Gf(p)). 

p a,6=± 

(22) 

In the limit of weak scattering 1/r 0, ImG''^ becomes 
strongly peaked at the Fermi energy. As a result intra- 
band contributions (a = b) dominate in this limit. Ne- 
glecting then the inter-band terms (a ^ b) and with the 
approximation 



ImGf ImGf « TTT5{e - Ea + fi) 



(23) 




one recovers the expression for the conductivity of the 
Boltzmann equation approach, when we identify Ttr with 
T. In the general case with a finite scattering rate, the 
inter-band terms and deviations of ImG^ImG^ from the 
delta function have to be taken into account. This is 
clearly demonstrated in Fig.|^. 
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FIG. 3. Act versus magnetic field for app = 0.4 and var- 
ious scattering strengths. Scattering events are taken into 
account in the broadening of the energy levels as discussed 
in the text. 

Taking into account a finite l/r the sharp edge of 
the anisotropy at Ug — app disappears, although the 
anisotropy curves do not change much even in the strong 
scattering case where epT = 0.8. 

Going beyond this crude approximation in calculating 
the conductivity we have to pay the price to make an 
assumption on the scattering mechanism at play. We 
assume elastic scattering, and we neglect quantum inter- 
ference corrections to the conductivity. In the following 
we will first calculate the conductivity in the case of a 
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pure s-wave scattering. The results will be obtained in 
the framework of the self-consistent Born approximation. 
Then we will allow an arbitrary p-dependence of the scat- 
tering potential. In that case we will confine ourselves to 
the weak disorder limit. 



IV. S-WAVE SCATTERING 

In this section we consider a specific microscopic origin 
for the conductivity, namely elastic scattering. For sim- 
plicity we assume a impurity potential with short range 
Gaussian correlations 



(L/(x)L/(x')) 



1 



27riVoT 



S{x - x') 



(24) 



We will calculate the conductivity in the framework of the 
self-consistent Born approximation, including the corre- 
sponding vertex corrections. 

The self-energy is then given by 



S = 



1 



2TrNnT 



■E 

p 



Gip). 



(25) 



To deal with the above matrix equation we expand it in 
Pauli matrices. The Green function can be written as 
G = Gqcto + GifTi + 6*2(72 with 



+ G_ 
apy 



Si 



apx 

2n 



2n 



G_ 



where il{p) — \J {apy 



a'^p^ and 



G± = (e - p^/2m + fiT^ - ^o)^^- 



(26) 



(27) 



Because the self-energy shares the matrix structure of the 
Green function, it has no (73 component. Furthermore it 
turns out that S2 = always solves the self-consistency 
equation, due to the symmetry p^ —Px- As a result the 
self-energy has the form Spcro -|- SiCTi. To appreciate the 
meaning of the two self-energy components, we briefly 
consider the case with no spin-orbit coupling. This is 
especially important for understanding the behavior at 
large magnetic field. The real part of the self energy Eg 
shifts the energy spectrum by a constant. Since we have 
to adjust the chemical potential in order to keep the par- 
ticle number fixed, ReSg can be safely neglected. The 
real part of Si may be re-absorbed into a renormaliza- 
tion of the Zeeman energy. We then concentrate on the 
imaginary part, and find 



ImS'^ = - 



1 

2^ 



■(70 



1 

2^ 



■(71 



(28) 



1 




To 


T 2No 


1 




n 


T 2No 



(29) 
(30) 



The sum or difference 1/Tip — 1 /tq ± I/ti are nothing but 
the scattering times for the two spin subbands. For weak 
magnetic field (wg < ef), the density of states in the two 
subbands are identical and therefore 1/t± = I/tq = 1/t. 
In the strong field limit (lus > ep) the upper band is 
depopulated, so that l/n = I/tq — l/2r. The scatter- 
ing rate in the empty band is zero due to the vanishing 
density of states. Typical numerical results are shown 
in Fig.^. Although the general shape of the curves is as 
discussed above, it is seen that for the disorder strengths 
we consider there are already considerable modifications. 
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4. The imaginary part of the self-energy for 
0.4 and epr — 0.2, 0.4 as a function of the mag- 
The numbers for epr = 0.4 have been divided by 



A. The vertex corrections 

We now move to consider vertex corrections, following 
the standard approach in the literature (see for exam- 
ple [^). Physically vertex corrections are important for 
obtaining the correct momentum relaxation time, which 
generally differs from the quasi-particle lifetime. For 
the choice of a pure s-wave scattering potential, as in 
Eq. (p4[) , in the absence of spin-orbit coupling, the mo- 
mentum relaxation time coincides with the quasi-particle 
lifetime. This corresponds, diagrammatically, to the fact 
that there is no dressing of the current vertex with mo- 
mentum independent impurity lines. In the presence of 
spin-orbit coupling this is no longer the case due to the 
presence of extra terms in the current vertex of Eq. (^7|) . 
The dressed current vertex is obtained by solving the 
equation 



jx,y _ ■x,y , 



1 



2TrNoT 



p 



(31) 



with 



A similar equation exists for the other energy ranges 
when one has a pair of retarded or advanced Green func- 
tions. The strategy to solve the equation is to sepa- 
rate momentum dependent from momentum independent 
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parts, expand the spin-dependent quantities in the Pauli 
matrix basis, and finally invert the resulting matrix equa- 
tions. 

We begin by separating the p-dependent and the p- 
independent parts of the current operator: 



f''" = Px,y/m =F aa2,i 



(32) 
(33) 



Since we consider a p-independent impurity scattering 
the Jq'^{p) part is not dressed, i.e.. 



(34) 



The equation for the momentum independent part of the 
current vertex reads 



ss' 'ss' 



1 



2TTNnT 



bs' 



(35) 



j:. 



2nNoT ^ 
p 



G«(p./m)G^ 



(42) 



One might try to determine the dressed vertices F^'^ 
by solving the above matrix equations. Whereas this is 
straightforward for Ff 3 , the equation for Fq needs spe- 
cial care. The analysis of Eq.(39) shows that the matrix 
on the right-hand side of (|3^) has an eigenvalue equal to 
one. As a consequence the matrix 

1 - (/oo + h\ + h2) -{loi + ho) 
-{lai + ha) 1 - (/oo + hi - ^22) 

cannot be inverted, and (^) has no unique solution. This 
property follows from charge conservation and manifests 
in the charge density vertex being singular in the low 
frequency, long wavelength limit (see also App.^). To 
see explicitly this property consider the imaginary part 
of the self-energy: 



We included here the spin indices s, s', a, and b 

The quantities 7^'^ are a sum of the bare vertices =Foo'2,i 

and a term which is generated by Px.y/m, 



Is 



Ta{cr2,i)s 



1 



J2J:—GfaGi,. (36) 



In the Pauli matrix space the vertex equation becomes: 



^7 + 



(37) 



2TrNaT ^ 



G^) 



p 



2'kNqt ^ 



G^(E«~I]-4)G'4. 



(43) 



By means of Eq.(^3|) we can relate the imaginary part of 
the self-energy with the integrals over pairs of retarded 
and advanced Green functions I^u- Expanding in Pauli 
matrices we find 



(E^ 



1 



S^)p = ^TT{aPa^'a^a'')I^,{ll'' - S-^)a (44) 



with 



which gives rise to the equation 



p 



(38) 



Some of the integrals /^i^ are zero and so the equations 
simplify. The final result is 



^ 2 
^ 3 



V 

7o 

V 

loa 
hi 



72 



ho 



hi 
ho 



I: 



22 hi 
ho 



ho 
hi 



111 + hi i(/oi ~ ho) 



^i(/oi ^ Iio) ho ^ 111 + ^22 

where the bare vertices 7^'^ are 

ryV — -fV t IV _L_ TV 

10 ^ '^00 ' "'11 ^ "^22 



7f = a- 

72^ 



•^01 



■^10 



-a - 



'02 



'20 



with 



Ff 



1/ro 
1/ri 



ho + hi + ^22 ^01 + ho 
hi + ho ho + hi ^ hi 



l/ro 
(45) 



This means that the vector of the scattering rates is an 
eigenvector of the above matrix with eigenvalue Aq = 1. 
The second eigenvalue is then given by Ai = 2/oo + 2hi — 
1 . We now demonstrate that the eigenvector correspond- 
(39) iiig to Ai is proportional to (7oi7i)- We start from 



— TrV/(G^-G'^) =0. 



2'kNqt 



(46) 



(40) 



We use again the relation G^ - G^ = G^\Yh^ -11^)G^ 
and arrive at 



1 



E 



4ttNot '■ — ' m 
p 



Tr 



1 

— ( 



Tl 

(47) 







(41) 



with the final result that 



To 



00 



Jli + ^12) + -(« + ^o'l + Jlo) = 0, (48) 
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so that the vector {-Jq, ^f) is perpendicular to (I/tq, l/n) 
and therefore necessarily is the second eigenvector of the 
matrix on the right hand side of (^5|). The solutions 
of ( ^5| ) for the dressed vertex {T^,T\) are then given 
by {T^,rl) = (7Q,7f)/(l — Ai) plus an arbitrary vec- 
tor proportional to (I/tq, 1/ti). From and ( p7| ) one 
observes that this arbitrary vector is irrelevant for the 
conductivity. For simplicity we choose it to be zero. 

Numerical results for the dressed vertices are shown 
in the figure. Notice that Eqs.(^) and ( |4^ ) provide a 
nontrivial consistency check for the numerics. 
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FIG. 5. The dressed vertices Fq . . . Fs in units of a as a 
function of the magnetic field. Here we chose app — OAep 
and l/(eFT") = 0.2. Remember that the bare vertices are 
7q = 73 = and = —7! = a; the asymptotic values in 
the strong field limit are Fj — — Fq = a/2, F2 — —a, and 
Fs = 0. Notice that F2 reaches the asymptotic value for 
high magnetic field only very slowly. 

It is seen that the anomalous part of the current op- 
erator is strongly affected by the vertex corrections. In 
particular for weak magnetic field the quantities F^ be- 
come extremely small, so that the anomalous part of the 
current operator is canceled, J* « p^/m. In the limits 
epT ^ 1 and lUs = this cancellation is exact as we now 
show. When tj, = the Green functions are 



G, 



R(A) 
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G 



B.{A) 



RiA) 



BiA) 



G 



R(A) 



2 

sin(6l) 
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G 

RiA) 
+ 

QRiA) 



G 



R(A) 



G 



R{A) 



(49) 



In the strong field limit oJs ^ 00 we find = 
Px/m - aa2, Jy = Py/m + a{ai -ao)/2, compare Fig||. 
Thus for Jx the vertex corrections are absent but they 
seem to be present for Jy. On the other hand as we 
noted before the solution of the vertex equation is not 
unique and also Jy ~ Py/m + aai solves (^), i.e. there 
are no vertex corrections to the conductivity. We will 
discuss the strong field limit in more detail in the next 
subsection. 



B. The conductivity 



Having discussed the vertex corrections we can now 
consider their effect on the conductivity. We do not find 
a strong magnetoresistance, compare also Fig.^ below. 
Typical numerical results for the anisotropy in the con- 
ductivity are shown in Fig.^. The most striking features 
are the change of sign at low magnetic fields and the 
sharp structure around lOs = £f- 
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FIG. 6. Anisotropy in the magneto-conductance for 
s-wave impurity scattering, l/{eFT) = 0.2 and various 
spin-orbit energies. For comparison we included the results 
of the relaxation time approximation {app = OAep', full 
line) . For not too strong fields, the vertex corrections change 
the sign of the effect. 

In Fig.0 it is seen that the structure around tOg ^ ep 
becomes more pronounced the less disordered the system 
becomes. For epT <^ 1 a step in the anisotropy evolves 
at that energy. The step is related to the van-Hove sin- 
gularity in the density of states at the band edge of the 
upper band. Disorder smears out this singularity. 



G 



RiA) _ 



± 



= (e - p /2m =F ap ± i/2T)" 



(50) 



The first step is to evaluate the effective bare vertices en- 
tering the expression for the conductivity, i.e., 7*. This 
can be done by evaluating the integrals we denoted by 
J^y. Because the Green function in the case oJs — 
has a simple angular dependence, one may directly con- 
clude that the integrals Jqo, Ju, and J22 are zero. With 
the approximation G^G^ ~ 2nTS{e — E±) the remain- 
ing integrals are readily found to be Jqi + Jio ~ ~a and 
J02 + J20 = ct with the result that the vertices 7* and 
consequently the dressed vertices F^ are zero. 
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FIG. 7. Anisotropy in the magneto-conductance for 
app — 0.4 and varying the disorder. The curve in the limit 
of weak disorder was obtained within the approximation de- 
scribed in the section on small angle scattering. 
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In order to have an understanding of these effects 
we will calculate the conductivity analytically for strong 
magnetic fields. The limit LOg —>■ oo is understood easily. 
There is just one band occupied with the dispersion given 
by E-{p) — /2m — Ws + ap^ . The Fermi surface is a 
circle which is centered around (0, —ma). Shifting p^ to 
p'v — pV fjiQ, problem becomes rotational symmet- 
ric. The (bare) current operator in the lower band is then 
simply of the standard form j^'^ — p'^-^ /m. As a con- 
sequence the p-independent disorder potential does not 
renormalize the current operator, and the conductivity is 
<7xx — o'yy = Co- particular there is no anisotropy in 
the conductance. 

Lowering the magnetic field various contributions to 
the anisotropy appear. Expanding the square roots in 
the energy dispersion, the Green functions for the two 
subbands read 



G^(c.,p) = 



where 



2to4 



Py ^ ^-1 



m± 



1 ± 



T± To Tl \ 2 ^2 



(51) 

(52) 
(53) 



We discussed the effective mass anisotropy already in the 
framework of the relaxation time approximation. Allow- 
ing the scattering times in the two bands to be different, 
the result of section || is generalized to 



1 (app)'^ T^ri- 



■ T+n+ 



2 uJs^F T^n^ -\- Tj^Uj^ 



(54) 



Further contributions to Act are due to the anisotropic 
relaxation rate. The scattering times tq and ti must be 
obtained by solving the self-consistent Born approxima- 
tion. In the weak disorder limit the van-Hove singularity 
at the band edge of the upper band leads to a singu- 
lar (=step) behavior of the two scattering rates. Such 
singular behavior is smeared out when disorder becomes 
stronger (compare Fig.^). For simplicity we restrict the 
following discussion to the weak disorder case. The an- 
gular dependent rate which we obtained by expanding 
the square root in the energy dispersion may also be ob- 
tained by considering the scattering probability in the 
eigenstate basis. The latter can be written as 



1 / 1 -t- 



pp 



1 ©p • Gp/ 



2 V 1 - ep ■ Gp/ 1 



(55) 



The vector Gp has been defined below Eq.(p!5|). In the 
high field limit 



off 



1 
1 



2I a^fe 



P'x? 



1 -1 
-1 1 



(56) 



i.e. the effective scattering probability is anisotropic. 
The anisotropy in the scattering rate is due to the term 
proportional to (apa;)^/^^ and is found as 5{\/t±) = 
T2i:\U\'^{apx/2u}sf{N+ - N_). For ojs < ep the two 
density of states are practically identical, so l/r± re- 
mains isotropic. Only for LUg > ep the scattering rate 
is anisotropic and is a source for the anisotropy in the 
conductivity, which we determine as 



'yy 



1 (qpf)^ 

4 uj^ ■ 



(57) 



We now consider the vertex corrections. To leading order 
in 1/lOs we need to calculate only the diagram with one 
impurity line. It is of the type 

^'^^^ = 1^ E ^■^(P)GpGp^PP'Gp'Gp'^'^(P')' (58) 

P:P' 

which one has to sum over the two bands. Because the 
scattering probability depends in the strong magnetic 
field limit only on the x component of the momentum, 
only axx is affected. We find by evaluating Eq. 



'vv 



CTO 



l{apFl_ 
2 uj'l 



(59) 



When ojs < ep the vertex correction cancels with the con- 
tribution due to the mass anisotropy since (n+ — n_)/n « 
LOsjtF^ so that the anisotropy vanishes. In the case 
LOa> tp the sum of the three terms we discussed gives 



1 {appf 3 {appf 
IXa/cro = h — 

I UJs€F 4 L02 



(60) 



in full agreement with the numerical results of Fig.(|6|). 

We conclude this section by noting that the order of 
magnitude of the anisotropy in the microscopic calcula- 
tion is the same as in the relaxation-time approxima- 
tion. Here however also the transport time becomes 
anisotropic. The competition with the mass anisotropy 
kills Act when Ws < cf- In the next section we will see 
that this cancellation is a peculiar consequence of the 
momentum independent impurity scattering. 



V. SMALL ANGLE SCATTERING 

We will now drop the assumption of a p-independent 
impurity potential and consider the more general case. 
For simplicity we will restrict ourselves to the weak dis- 
order limit, where we will neglect the disorder broadening 
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of the energy levels, and terms which are related to the 
off-diagonal matrix elements of the current operator (in 
the eigenstate basis). The results obtained within these 
approximations are fully equivalent with the Boltzmann 
equation approach, as it is discussed in Appendix B. In 
the energy state eigenbasis the dressing equation for the 
current vertex then reads 

jx,y _ -x.y I Ti^eff r<R r<A rx,y z^-, 

pni Jpm ^ / J ' pm-jp' m' ^ p' m' ^ p' 7n' ^ p' 7n' \^ / 
p'ra' 



with 



(62) 

The lifetime is determined from the expression 



We determine the dressed current operator with an ex- 
pansion in eigenfunctions, following Ref. [Q: let (/)p„ be 
the full set of eigenfunctions of the operator 



ETi^cff f^R 1 1 _ J. 1,1 

pm\p' ra' p' m' p' ra' T p' m' V^pm ' 

p'ra' 

which are normalized according to 

pm 



(64) 



(65) 



The kernel entering the integral equation ( pl|) may be 
expanded as 



(66) 



One further notices that l/rpm is an eigenfunction with 
eigenvalue one, A*^ = 1 and oc 1/Tpm. Furthermore 
the bare current operator is perpendicular to this eigen- 
function. This can be shown in analogy with the discus- 
sion developed in the case of s-wave scattering in section 
IV. A. With the expansion 



3 



pm~/_^Jl Vprrn Jl ~ /_^Jpm'^pm^pm'rpm- V" ' 7 



/ pm 

the dressed current operator can be expressed as 



J. 



pm — 2^ Vpm 1 _ \ 



3i 



x,y 



(68) 



from which wc finally obtain the conductivity by insert- 
ing it in Eq.(j2^). Before showing the results for the 
anisotropy and analyzing the vertex corrections in the 
presence of small angle scattering, it is instructive to 
consider the zero magnetic field case first. As discussed 
in section IV for the case of s-wave scattering, the ver- 
tex corrections have their most dramatic effect on the 



anisotropy in the limit of weak magnetic field. In the zero 
magnetic field limit, the structure of the dressed current 
vertex simplifies considerably. The scattering probabil- 
ity Wpp' and the effective scattering probability become 
isotropic, i.e. they depend only on the angle between p 
and p'. One can then expand 



yycff ^ ^0 _^ 214^1 cos(6' - 9') + 



(69) 



where W^^, W^, and have still the two by two matrix 
structure. The lifetimes and the dressed current operator 
are expressed as 



1/Tp™ = 2nJ2 Km;p'm'6{f^ " ^^^'(p'))- (63) ^ud 



jx_,y 



1/t+ = 2TrN+Wl+ + 2TrN^Wl_ 
1/t_ = 2nN+W"_+ + 2ttN^W\ 



2ttN+t+WI^ 2ttN^t^WI^ 



(70) 
(71) 



2ttN+t+W}_+ 2ttN^t^W^_ J \ J- 



(72) 



where -/V±(e) = Ao(l T ma/ {ma)^ + 2me) and j± = 
Vp£'±(p) are the density of states and the bare current 
vertices in the two bands. For short range impurity po- 
tential, from Eq.(|5|), one finds W^^ = ... ^ I/AttNot, 
T± = r, and iy|+ = -T4^|_ = . . . = I/SttNot. Go- 
ing through the algebra one may verify our earlier re- 
sult that the vertex corrections cancel the anomalous 
velocity operator in this limit, 3± = p±/m. Gener- 
ally however this is not the case. In the extreme for- 
ward scattering limit in particular, one can neglect the 
off-diagonal terms and arrives at the standard ex- 

pression for single, isotropic bands J± = j±Ttr/T with 
l/rtr = N± J deWlli0)[l - cos(0)]. 

Let us now consider a Gaussian scattering probability. 



I^pp, - W^exp(-r2(p-p')V2). 



(73) 



The electron spin is assumed to be conserved, so this scat- 
tering potential has still to be transformed to the energy 
eigenstate basis in order to obtain W^p, , see Eq. ( |55| ) . The 
parameter rg controls the decay of the scattering poten- 
tial, tq = corresponds to pure s-wave scattering, and 
with toPf 3> 1 we are in the strong forward scattering 
limit. For simplicity the strength of the scattering po- 
tential and the screening parameter tq are here assumed 
to be independent of LUg. Of course, in order to explain 
the large magnetoresistance such a dependence has to be 
taken into account. Suggestions for a microscopic origin 
of the magnetoresistance have been given, e.g., in Refs. 
pl|-p^ . At the end of this section we will also consider 
the charged impurity scattering, which was discussed in 
the absence of the spin-orbit coupling, e.g., in p2| , ^ . 

Although our main concern here is the anisotropy in 
the magnetoconductance, it may be useful to begin by 
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showing the numerical resuhs for the magnetoconduc- 
tance itself (Fig.||) for the model with a Gaussian scat- 
tering probability. 



overall amplitudes of the various curves are scaled in or- 
der to make them more comparable. Also as a function 
of spin-orbit energy the position of the peak changes only 
weakly. 
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FIG. 8. Magneto-conductance for app — 0.3 and two dif- 
ferent screening parameters ro. Within the approximations 
in this section the conductivity is proportional to the inverse 
of the scattering strength, W, which we assume here to be 
independent of the magnetic field. For s-wave scattering 
(ro = 0) the anisotropy of the conductivity is not seen on 
the scale of the plot. For ropp = 1 we show the conductivi- 
ties parallel and perpendicular to the magnetic field. 

Given a scattering strength W, the conductivity in- 
creases with the parameter rg and with the magnetic 
field. Both effects are understood easily. The kine- 
matics selects scattering processes with q = |p — p'| ~ 
2pF, whereas the scattering potential of Eq.j?^) gives 
q < I/vq. The condition for effective scattering becomes 
ToPf ^ 1- Increasing 7'o weakens the above condition 
and the conductivity gets enhanced. As a function of the 
magnetic field the conductivity increases due to a similar 
reason. The Fermi surface of the lower band grows and 
the effective condition ropp.- ^ 1 is no longer satished. 

The numerically determined anisotropy in the magne- 
toconductance is shown in Fig.^. Since now the conduc- 
tivity itself is strongly magnetic field dependent we scale 
Act with axx + o'yy instead of ctq. 
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FIG. 9. Anisotropy in the magneto-conductance 
app = 0.3ep and varying the screening parameter 
The screening parameters ropp = 1,2,4 correspond 
Ttr/r ~ 2, 8, 30 when they are determined for app — 



for 

ro- 
to 
= 0. 



One observes that the p-dependent scattering increases 
the overall size of the anisotropy up to a factor ~ 5 com- 
pared to the s-wave case. The sign of the anisotropy also 
changes. The peak anisotropy is still at magnetic helds 
of the order of the Fermi energy, but there is a weak de- 
pendence on the scattering potential. The anisotropy as 
a function of spin-orbit energy is shown in Fig.nn. The 



-0.2 
-0.25 




FIG. 10. Anisotropy in 
strong forward scattering 



the magneto-conductance for 
{ropp = 4) for different 
strengths of spin-orbit scattering. The curves for 
app/ep = 0.2, 0.1, 0.05 are multiplied by 4, 8, 16. 



Finally let us consider charged impurity scattering. 
This type of scatterin g m ay explain certain aspects of 
the magnetoresistancc |32|] . For a charged impurity situ- 
ated in the plane the scattering potential is 



U{q) = 27reV[ge(g)], 



(74) 



where e{q) is the dielectric function. We approximate it 
here as 



1 



1 



1 + (2^e2/g)x(9) 



(75) 



where xil) is the two-dimensional Lindhard function It 
has been argued in Ref. |^ that the screening wave num- 
ber Qs = 'iire'^ N{eF) is for low electron density large com- 
pared to the Fermi wavelength with the consequence that 



Uiq) 



xiq)' 



(76) 



The scattering probability is approximated according to 
W{q) cx \Uiq)\\ 

In Fig.|ll| we show the resistivity obtained with app = 
0.2 and 0.4eF- 
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FIG. 11. Resistivity versus in plane magnetic field for 
charged impurity scattering. 
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The resistivity increases by a factor of four from the 
zero field hmit to very strong magnetic fields. In the ab- 
sence of spin-orbit coupUng the resistivity saturates when 
uis > ef, where only one Fermi surface is left |3^. In the 
presence of spin-orbit coupling there remains a weak in- 
crease of the resistivity up to higher magnetic fields, since 
the depopulation of one band is no more equivalent with 
the full spin polarization. 

The anisotropy in the magnetoresistance for charged 
impurity scattering is shown in Fig.O. 



0.05 
0.04 
0.03 
0.02 
0.01 


-0.01 
-0.02 
-0.03 
-0.04 





1 


1 


- apF 


= 0.2ejj 




- apF 


= 0.4eF 






1 1 1 



0.5 



1 



FIG. 12. Anisotropy in the magnetoconductivity versus 
in plane magnetic field for charged impurity scattering. The 
strength of the spin-orbit field is app = 0.2, 0.4eF. 

Within the approximation adopted the scattering 
probability depends only weakly on the momentum. 
Therefore results for the anisotropy are in the interme- 
diate way between the s-wave and the strong forward 
scattering limit. 



VI. DISCUSSION 

We have calculated the Drude conductivity for a two- 
dimensional Fermi gas with spin-orbit coupling. In the 
presence of an in-plane magnetic field the conductivity 
parallel and perpendicular to the field are different. 

Such an anisotropy has recently been observed in Si- 
MOSFETs m. For the range of densities in the ex- 
periment Us w 0.7 — 1.3 X 10^^/cm^, the Fermi energy 
is ep ~ 5.1 — 9.5K. Taking a « 6 x 10~^Kcm which is 
reported in Ref . ||2^ , the dimensionless spin-orbit param- 
eter is of the order app/^F ~ 0.3 — 0.7. From our theory 
one expects a relative anisotropy in the magnetoconduc- 
tance of the order of several percent, in agreement to 
what is seen experimentally. The maximum anisotropy 
appears for magnetic fields of the order tOs ^ ep, so the 
peak energy scales with the density, which again agrees 
with the experiments. We found that the details of the 
effect depend on the scattering potential. In the case 
of pure s-wave scattering the maximum anisotropy ap- 
pears precisely at the magnetic field, where the upper 
band depopulates. The conductivity is larger when it is 
measured parallel to the magnetic field. On the other 
hand, at higher fields, lUs > ep, the conductivity is larger 
when measured perpendicular to the magnetic field. The 
sign of the anisotropy at low fields is at odds with the 



experiment. Allowing a finite range of the scattering po- 
tential, we have shown that the anisotropy sign changes. 
In particular, in the case of small-angle scattering both 
the sign and the peak position agree with what is found 
experimentally. For the case of charged-impurity scatter- 
ing we found - within the random-phase-approximation 
- a rather complex pattern of the anisotropy. There are 
two peaks with opposite sign, one peak near tUs ^ app 
and the other at lJs ~ ep- These results however should 
be taken with caution when comparing with experiments 
for the low density electron gas. 

Comparison with experiments in heterostructures is 
even more delicate due the intrinsic crystallographic 
anisotropics. In both the experiments reported in ^5 26 
the anisotropy was measured along two crystallographic 
directions with different mobilities and most likely with 
different operating scattering mechanisms. At low den- 
sities (close and across the MIT) the anisotropy shows 
Act < 0, while at high densities there is a change of sign 
as function of the magnetic field: Act < (>)0 for low 
(high) fields. In our opinion, while a comparison with 
the sign of the anisotropy contains the uncertainty re- 
lated with the knowledge of the scattering mechanism, 
the fact that the maximum anisotropy at low and mod- 
erate fields scales with the density appears as a robust 
feature in agreement with our theory. 

In summary we think that the physics described in this 
paper may explain the anisotropic magnetoresistance ob- 
served in Ref. in Si-MOSFETs, and perhaps that in 
heterostructures |25| , |2^ ■ At more general level, our re- 
sults give a further hint that the Drude-Boltzmann the- 
ory is a good starting point for the description of the 
transport properties of the metallic phase of the 2D elec- 
tron gas, even near the observed transition to an insula- 
tor. 

We acknowledge fruitful discussions with C. Castellani 
and V. Falko. This work was supported by the DFG 
through SFB 484 and by EU Research Training Network 
program (Project Nr: RTNl - 1999-00406). 



APPENDIX A: WARD IDENTITY 

In this appendix we derive for the sake of complete- 
ness the Ward identities for the system with spin-orbit 
coupling. The generalized continuity equation is 



atp(x,i) + a,r(x,i) = 



(Al) 



where i = x,y, z. The density and the current operators 
are 



p(x,t) = ^Z-t (x,t)'(/'a(x,i), 



r(x,t) = ^[(9,V'i(x,i))^a(x,i) 



(A2) 
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V'i(x,t)(9,Va(x,t))] 

aei-jk-rijipl (x, t)ak.af3'ipt3{'^, t). 



(A3) 



Define now the vertex functions 



A°^ = (rip'(x,t)V'a(x',OV'^(x",i")) 
A^„^ = (rir(x,i)v^„(x',t')^^(x",<")), 

for wliicli tlie Ward identities read 

= 6{t - t")6{x - x")iG„/3(x't'; x, t) 
- d{t - t')6{yi - x')iGa,3(xt; x"t"). 

After Fourier transform one gets 

In the static limit lo = Q and letting q one gets 

d 



(A4) 



(A5) 



(A6) 



= -Gc,a{p,e)Jla'{P:0)Ga.p{p,e) 



dp' 



(A7) 



The above equation is sufficient to see the cancellation of 
the diamagnetic term in the conductivity, see Eq.(n9|). 



^Tr[j\p)G{p)J^{p)G{p)]^^5., 



— oo 
oo 



^E / Ott'^'^ 



27r 



-Tr 



np)-^Mp) 
^.mG{p) 



N 



N 



— 5v 



1 r°° <\f 



N 



(A8) 



In the dynamic limit on the other hand g = 0, a; — > 0, 
Eq.(O) becomes 



APPENDIX B: BOLTZMANN EQUATION 

We now briefly demonstrate the equivalence of the ap- 
proach in section ^ with the Boltzmann equation. Since 
we have in mind to solve the Boltzmann equation in the 
presence of a DC electric field, the distribution function 
gpa is chosen to depend only on the wave vector p. The 
Boltzmann equation including elastic scattering is then 
written in the form 

- eE • Vgpa = - ^ (9pa,p'/3 (ffpQ - gp'p) ■ (Bl) 
p'/3 

The Greek indices a and /? label the two spin subbands. 
The scattering kernel, in the case of elastic scattering, is 
related to the scattering probability by 

Qpa,p'/3 = 2ttS (i?p„ - i?p,^) W;l^,f, (B2) 

A solution of the linearized Boltzmann equation is looked 
for in the form 



df 

gpa = f{Epa) + ■ Upa, 



(B3) 



where the vector function Up^ obeys the integral equa- 
tion 



Vpa = ^ Qpa,p't3 (Upa — Up'^) . 

p',/3 



(B4) 



Finally, the electrical current density is given by 

j = -e^Vpa5pa, (B5) 



which becomes 



j = ^'E(-|^)^p"("p"-E)- 



pa 



(B6) 



To make contact with the diagrammatic approach we 
write 



(A9) 



which can only be solved if the density vertex A'^ is singu- 
lar in the zero frequency limit. To see how this translates 
for the irreducible vertex F, write A°^ = Gaa^ aa'Ga'p 
and perform the p-summation on both sides of Eq.(A9). 
After decomposing in the Pauli matrices components, one 
gets 



ro 
ri 



70 
71 



l/ro 



(AlO) 



where we have used Eqs.(|39D. From Eq.( [AlC| ) one sees 
that both Fq and Fi have to be singular in the zero fre- 
quency limit. 



with 



Upa — Jpa'^po 



p'l3 



(B7) 



(B8) 



The vector function Jp^ is the renormalized current ver- 
tex, which arises in the diagrammatic approach. This is 
easily seen in the T = limit, by observing that 



^ ^ - Epa) ~ T^;;;— Gp„Gp„ 



27r T-pa 



(B9) 



Upon using Upo, = JpaTpc and Eq.(B9), the integral 
equation (B4) becomes 



12 



(BIO) [15 



Hence the linearized Boltzmann equation is equivalent 
to the approach in section where we neglect lifetime 
broadening. 
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